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Abstract 



Using algebraic Bethe ansatz and the solution of the quantum inverse scatter- 
ing problem, we compute compact representations of the spin-spin correlation 
functions of the XXZ-^ Heisenberg chain in a magnetic field. At lattice dis- 
tance m, they are typically given as the sum of m terms. Each term n of this 
sum, n = 1,...TO, is represented in the thermodynamic limit as a multiple 
integral of order 2n + 1 ; the integrand depends on the distance as the power m 
of some simple function. The root of these results is the derivation of a com- 
pact formula for the multiple action on a general quantum state of the chain 
of transfer matrix operators for arbitrary values of their spectral parameters. 
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1 Introduction 



The main challenging problem in the field of quantum integrable models is to compute exact 
and manageable expressions for their correlation functions. This issue is of great importance, 
not only from a theoretical view point but also for applications to relevant physical situations. 

The archetype of quantum integrable lattice models is provided by the XX Z spin-^ Heisen- 
berg chain in a magnetic field, 

M 

H=Y. + «+i + - 1)) - hS„ (1.1) 

m=l 

where 

1 

•5- = 2E^-' [H,S,]=0. (1.2) 

m=l 

Here A is the anisotropy parameter, h an external classical magnetic field, and a^'^ denote 
the usual Pauli matrices acting on the quantum space at site m of the chain. The simultaneous 
reversal of all spins is equivalent to the change of the sign of the magnetic field, therefore it is 
enough to consider the case /i > 0. In the thermodynamic limit M ^ oo and at zero magnetic 
field, the model exhibits three different regimes depending on the value of A: for A < —1, the 
model is ferromagnetic; for — 1 < A < 1, the model has a non degenerated anti- ferromagnetic 
ground state, and the spectrum is gapless (massless regime); for A > 1, the ground state is 
twice degenerated with a gap in the spectrum (massive regime). 

Although the method to compute eigenstates and energy levels goes back to H. Bethe in 1931 
H, m, ^, the knowledge of its spin correlation functions has been for a long time restricted 
to the free fermion point A = 0, a case for which nevertheless tremendous works have been 
necessary to obtain full answers ||5|, |6|, 0, 10 1. 



However, at zero temperature and for zero magnetic field, multiple integral representations of 
elementary blocks of the correlation functions (see definition below) have been obtained from the 
g-vertex operator approach (inspired from the corner transfer matrix technique) in the massive 
regime A > 1 in 1992 |]l^], and conjectured in 1996 |]l^] for the massless regime —1 < A < 1 
(see also |13]). A proof of these results together with their extension to non-zero magnetic field 



has been obtained in 1999 ^ for both regimes using algebraic Bethe ansatz ^] and 
the actual resolution of the so-called quantum inverse scattering problem [0, |l^ . 
These elementary blocks for correlation functions are defined in the following way: 
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Here \ipg) denotes the ground state in the massless regime or any of the two ground states con- 
structed by algebraic Bethe ansatz in the massive regime; E^'^™' are the elementary operators 
acting on the quantum space Hm at site m as the 2x2 matrices E^j^ = 81^^18^^^. Any n-point 
correlation function can be reconstructed as a sum of such elementary blocks. 

To compute these elementary blocks, the following successive problems have to be addressed 
Ip^ , p^ : (i) determination of the ground state {ipg\, (ii) evaluation of the action of the product 

of the local operators E^^' ^ on this ground state, and (iii) computation of the scalar product 
of the resulting state with 

For the XX Z spin-i Heisenberg chain in a magnetic field, these problems have been solved 
in [0, 15| in the framework of algebraic Bethe ansatz. The central object of this method is 
the so called quantum monodromy matrix depending on a complex variable A (the spectral 
parameter). For the XX Z spin-^ chain, it is a 2 x 2 matrix with operator valued entries acting 
in the quantum space of states H: 

The quadratic commutation relations between these four operators are given by the Yang-Baxter 
algebra. It is governed by a trigonometric R matrix solving the Yang-Baxter equation. The 
Hamiltonian of the chain is then contained in the commutative family of operators generated 
by the transfer matrix T(A) = (A-\- D){X) for arbitrary values of A. The algebraic Bethe ansatz 
leads to the simultaneous diagonalization of these transfer matrices and of the Hamiltonian. 
The ground state {ipg\ (resp. \ipg)), as the other eigenstates, is given as the successive action 
of operators C(Afc) (resp. i?(Afc)) on the ferromagnetic reference state (0| (resp. |0)) with all 
spins up. Namely, we have {^pg\ = {0\Y\f^C{\k) and \Tpg) = Y\i^B{Xk)\0) for a particular set of 
spectral parameters {A^} solving the Bethe equations. 

To evaluate the action of local operators on this state, the strategy is to imbed them in 
the Yang-Baxter algebra of T matrices by solving the quantum inverse scattering problem (see 



H] for details) as 

Ef'^' = Ui^ + m)Te,4V fliA + Dr'Cj), (1.5) 

k=l k=l 

where coshr? = A. Then, using the Yang-Baxter algebra, one can reduce any elementary blocks 
of the correlation functions to multiple sums of scalar products of some states with {i/jg}- Each 



of these scalar products can be computed as the ratio of two explicit determinants 14]. In 



the thermodynamic limit, these multiple sums lead to m fold integrals over contours Cj which 
depend on the value of j, on the regime considered and also on the value of the magnetic field. 



The answer can be written generically as |15] 

m „ 

Fm{{e„e'^}) = n / dX, 17„({A}, {e„e'^}) 5,({A}). (1.6) 
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Here r2m({A}, {ej,e'j}) is a purely algebraic quantity, which in particular does not depend on 
the regime nor on the magnetic field. In contrast, 5/i({A}) is a functional of the density function 
/0/i(A) solution of the Lieb equation describing the ground state Q, and hence depends both on 
the regime and on the value of the magnetic field h. 

It is remarkable that, for zero magnetic field, the two multiple integral representations 
resulting from the q-veitex operator approach and from Bethe ansatz are identical: contours 
and integrands coincide. It would be very desirable to understand this intriguing fact directly 
at the operator level. Note that for non zero magnetic field, the quantum affine symmetry used 
in the (7- vertex operator approach is broken, and no result is known up to now from this method 
in this case. 

In principle, any n-point correlation function can be obtained from these elementary building 
blocks. One should note however that, although these formulas are quite explicit, the actual 
analytic computation of these multiple integrals is missing up to now. Moreover, the evaluation 
of correlation functions of physical relevance, as for example the spin-spin correlation functions 
at distance m on the lattice like {afa^^i), is a priori quite involved. Indeed, the identity 

m 

shows that the corresponding spin-spin correlation function is actually given as a sum of 2*""^ 
elementary blocks. So, the number of terms to sum up grows exponentially with m, making the 
problem of asymptotic behavior at large distance extremely difficult to solve in these settings 
from the present knowledge of the elementary blocks (|l.6| ). In the language of algebraic Bethe 
ansatz, and using the solution of the quantum inverse scattering problem, this question amounts 
to the computation of the following average value: 

(V^,|C7(2)(^ + Z)r-l(2)i?(2)|V,,). (1.8) 

Hence, to obtain manageable (re-summed) formulas for spin-spin correlation functions, and to 
avoid the computation of the above sum of 2"^"^ terms, we need to derive a compact expression 
for the action of the shift operator {A + D)"^~^{^) (from site 1 to site m + 1) on arbitrary states. 

The main purpose of this paper is to give a solution to this problem in the framework of 
algebraic Bethe ansatz, and to apply it to the evaluation of spin-spin correlation functions. 

In fact, for later use, we will solve an even more general question: the evaluation of a compact 
formula for the multiple action of transfer matrix operators {A + D){xa), for any set of spectral 
parameters Xa, on arbitrary quantum states (a priori not eigenstates) of the XX Z model (see 
Proposition |4.1| ). This leads to the evaluation of the spin-spin correlation functions at lattice 
distance m as the sum of only m terms (instead of 2"^~^), the n^^ term in the sum being 
expressed in the thermodynamic limit as a multiple integral of order 2n -|- 1, for n = 1, . . . , m 
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(see Proposition 6.1, Proposition |6.2|). For the two point correlation functions, a typical form 



of the result is ( |6.13D , 

K^^+i) = E fd-^'z I rX I [/({A,z}r ({A,^,z}) 5,({A,z}). (1.9) 

Here, a,/3 = x,y,z, and the functions F^'^dA, ^, z}), /({A,2;}) are purely algebraic quantities, 
which do not depend on the regime nor on the magnetic field; the integration contours C^, Ca, 
and the functional Sh of the density function describing the ground state (evaluated at points 
A and z) depend both on the regime and on the value of the magnetic field. 

In this formula, one can interpret the integrals over the A and z variables to be generically 
associated to the shift operator from site 1 to site m + 1, while the integrals over the variables 
fi correspond to the contribution of the operators cr",(T^. 

Hence this method provides us with an effective re-summation of the previous 2"^"^ elemen- 
tary blocks, although we found it more convenient and general to work it out at the operator 
level and in the algebraic Bethe ansatz framework. In particular, all our considerations are valid 
for the finite lattice case. We also believe that it can be applied to many other models for which 
multiple integral representations of elementary blocks of correlation functions are known, like 
for example the integrable higher spin Heisenberg chains [21|. 



Moreover, it should be stressed here that, for each term n of this sum, the distance appears 
now explicitly in the integrand merely as the power m of some function /({A,z}) of the inte- 
gration variables. This feature, which is the result of our re-summation, is obviously of great 
importance for future asymptotic analysis at large m. Let us finally mention that these new 
representations of the spin-spin correlation functions (valid for arbitrary values of A > — 1) 
lead in a simple way to the known results at the free fermion point A = 0; the corresponding 
computations will be presented in a separate publication. 

This article is organized as follows. In the next section, we recall some basics about the 
study of the XXZ spin-^ Heisenberg chain in the algebraic Bethe ansatz framework. In section 
^, we present the list of formulas necessary for the computation of the correlation functions via 
algebraic Bethe ansatz. In section we derive a compact formula for the multiple action of the 
transfer matrix operator on an arbitrary quantum state of the chain and for any value of the 
spectral parameters. It leads in section ^ to the evaluation of the generating functional of the 

cQj-j-giation functions. General spin-spin correlation function at lattice distance m are given 
in section |6|. Some perspectives are discussed in the conclusion. Lengthy computations and/or 
proofs of intermediate results are presented in a set of three appendices. 

We dedicate this paper to the memory of our friend and colleague A. Izergin. When we 
began this work two years ago, he was about to join us, but unfortunately these plans were 
suddenly stopped. 
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2 The XXZ spin-i Heisenberg chain 



The Hamiltonian of the cyclic XXZ chain with M sites is given by In the framework of 

algebraic Bethe ansatz, it can be obtained from the monodromy matrix T{X), which is in turn 
completely defined by the i?-matrix. The i?-matrix of the XXZ chain acts in the space (8) 
and is equal to 



RiX) 



1 



sinh(A + rj) 



( sinh(A + r?) \ 

sinh A sinh r] 

sinhr/ sinh A 

\ sinh(A + 7/) / 



coshr? = A. (2.1) 



It is a trigonometric solution of the Yang-Baxter equation. Identifying one of the two vector 
spaces of the i?-matrix with the quantum space Tirm one defines the quantum L-operator at 
site m by 

^m(A) = i?om(A-r//2). (2.2) 

Here i?om acts in ® 'Hm- The monodromy matrix T(A) is then constructed as an ordered 
product of the L-operators with respect to all the sites of the chain: 



r(A) 



^(A) B{\) 
C{\) D{\) 



Lm(A)...L2(A)Li(A). 



(2.3) 



The Hamiltonian (|l.l| ) at /i = can be obtained from T(A) by the trace identity 



H = 2sinhr/ — logr(A) 



const. 



Here, the transfer matrix 



T(A) =trr(A) = A{X) + D{\) 



(2.4) 



(2.5) 



generates a continuous set of commuting conserved quantities. For technical reasons it is con- 
venient to consider the inhomogeneous XXZ model, where 

Lm{\) = Lm{\U) = Ihm{\-U), T{X) = Lm{X,^m) ■ ■ ■ L2{X,a)Li{X,^i), (2.6) 

and are arbitrary complex numbers attached to each lattice site that are called inhomogeneity 
parameters. In the homogeneous limit S,m = ''?/2, we come back to the original model (|l.lD . 
The commutation relations between the entries of the monodromy matrix are given by the 
Yang-Baxter quadratic relation, 



RuiXi - A2)Ti(Ai)r2(A2) = T2(A2)Ti(Ai)fli2(Ai - A2). 



(2.7) 



The equation ( ^?7| ) holds in the space Vi (8) V2 (where Vj ~ C^). The matrix Tj{X) acts in 
a nontrivial way in the space Vj ®7i, while the i2-matrix R12 is nontrivial in Vi ®V2- 
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The space of states is generated by the action of creation operators B{X) and annihilation 
operators C(A) on the reference state |0) with ah spins up. In the following, we will consider 
general states of the form 

N 

\i;) = l[B{Xj)\0), N = 0,1,...,M, (2.8) 

i=i 

which are eigenstates of the transfer matrix T(fi) (and thus of the Hamiltonian in the homoge- 
neous case) when the parameters Xj satisfy the system of Bethe equations 

smh(Aj - Cm + 7?) fj^ smh(Aj - Afc - r/) 

The corresponding eigenvalue r(/i, {Aj}) of the operator T{fi) (containing the energy level fol- 
lowing (^)) is 

N N 

t(/., {A,}) = a(/.) J] /(A„ /i) + difi) n fif^, A,). (2.10) 
i=i i=i 

Here and further we use abbreviated notations for certain combinations of hyperbolic functions: 

ff. ^ sinh(A - /i + r?) sinh?? 

/ A,/x = — — TT hA,// = — — . . ■ ^ — r. (2.11) 

smh(A — fi) smh(A — fi) smh(A — fi + rj) 

The functions d{^) and a(/i) are eigenvalues of the operators D{^) and A{^) on the reference 
state: 

M 



d{^i)=\{r\^^,U), a(^) = l. (2.12) 



m=l 



Following the paper |15], we consider below the action of the monodromy matrix elements on 
the dual state which can be constructed similarly to (|2.8D via the operators C(A) as 



N 

(V| = (0|nC(Ai), iV = 0,l,...,M. (2.13) 

i=i 

Here (0| = lO)"*", and ( p. 131 ) defines a dual eigenstate if the parameters Aj satisfy the same 
system of Bethe equations 



Our final goal is to compute the correlation functions in the ground state in the thermody- 
namic limit M ^ oo. The thermodynamics of the XX Z chain was studied in j^, ^ Here 
we merely recall the formulas we need for our study. 

The ground state \ipg) of the infinite chain can be constructed as the limit of the finite chain 
eigenstate Uj=i B{Xj)\0) for M ^ oo, TV ^ cx) and N/M equal to some constant whose value 
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depends on the magnetic field h. In this limit, the Bethe equations for the set of parameters 
{Aj} reduce to the integral Lieb equation for the ground state spectral density ptotW'- 



c 



where 



sinh 27] 



(2.14) 



(2.15) 



sinh(A + rj) sinh(A — rj) 

The integration contour C = [— A/j,A/j] in ( |2.14 ) depends on the regime considered. In the 
massless case — 1 < A < 1, the contour C is an interval of the real axis and the parameter r] is 
imaginary: r] = —iC, C > 0. In particular, at /i — > 0, A/^ — > oo, and the Lieb equation can be 
solved explicitly (see ( p.l7| )). For A > 1 (r/ < 0) the limits =bA/j are imaginary, which means 
that the integral in ( |2.14| ) is taken over an interval of the imaginary axis. At h = 0, = —m 12 
and the solution of the Lieb equation is given in terms of theta- functions (see (|2.17| )). 

For technical purposes, we also introduce the inhomogeneous density /9(A, ^) as the solution 
of the integral equation 



2^ip{\i) 



[ K{\ 

Jc 



(2.16) 



It coincides with ptot{^) at ^ = ri/2. For our goals it is enough to consider —Q < Ini(^) < for 
— 1 < A < 1 and r] < Re(^) < for A > 1. For zero magnetic field, one has 



p(A,e) 



2Csinh|(A-0 



n 

n=l 



^2n 



1 + 



^2n 



^3(^(A-0,9) 



|A| < 1, c = in, 



A>1, 



(2.17) 



In the presence of the magnetic field, the equation ( p.l6| ) cannot be solved explicitely in terms 
of known elementary or special functions. Some properties of p(A, ^) can nevertheless be es- 
tablished : in particular, it is not difficult to see that p{\,i) has a simple pole at A = with 
the residue 27ri Res p(A, ^)|A=g = —1- This property was used in |15| for the computation of the 
elementary blocks, and we shall also use it in Section |6|. 



3 From quantum inverse scattering problem to correlation func- 
tions via algebraic Bethe ansatz 



In this section, we review the main steps of the method proposed in |14, 15] for the computation 
of the correlation functions in the framework of algebraic Bethe ansatz using the solution of 
the quantum inverse scattering problem: we first recall briefly how to obtain a multiple integral 
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representation for the elementary building blocks (O) (see [jl^, for details), then we discuss 
in this context the case of the spin-spin correlation functions at lattice distance m. 



The elementary blocks of the correlation functions are defined as the normalized expectation 
values of products of local matrices E^^' ^ from site j = 1 to site j = m with respect to some 
eigenstate of the transfer matrix (for which the set of spectral parameters satisfy the Bethe 
equations) : 

N 

(oi n c(Afc) 



k=l 



n 



N 

n BiXkM 

k=l 



N N 

(oi n c{\k) n BiXkM 

k=l k=l 



(3.1) 



For technical reasons, it is convenient to achieve the calculation in the generic inhomogeneous 
case ( |2.6| ); it is easy at the end to particularize the result to the homogeneous chain ( |1.1| ). 
To compute the expectation values (3.1), or more generally any kind of correlation functions, 

one has first to express the elementary local operators E^' ^ (or equivalently the local spin 
operators) in terms of the entries of the quantum monodromy matrix. Such a representation is 
given by the solution of the quantum inverse scattering problem |14, p!^ ]: 



Theorem 3.1. [14-, li] Let us consider the inhomogeneous XXZ model (2.t) with arbitrary 
inhomogeneity parameters S,k, ^ ^ k < M. The local spin operators at any site j of the chain 
can be expressed in terms of the elements of the quantum monodromy matrix as 



llT{Ck)-B{C,).llT-\Ck), 



k=l 



i-1 



k=l 
j 



n m) • c(Ci) • n ^^"'(6 



(3.2) 



k=l 



i-1 



k=l 



k=l 



k=l 



In particular, these formulas apply for the homogeneous model ^l.l\ ) where = ??/2, 1 < A; < M. 



Remark 3.1. The identity operator in the site j can also be written in a form similar to ( p.2D : 



\[T{i,).{A + D){i,).\{T-\i,). 



(3.3) 



fc=i 



fc=i 



Remark 3.2. The resolution (3.2)-( p73[) of the quantum inverse scattering problem can also be 

expressed in terms of the elementary matrices E-^' ^ in the site j. In the homogeneous case, the 
corresponding reconstruction formulas are given by (|1.l 
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From Theorem 3.1 and Remark 3.2, one obtains 



Fm{{e,,e'A) = ^^{{\])- 



(3.4) 



where $m({A}) is the ground state eigenvalue of the corresponding product of the transfer 
matrices: 



m N 



'fm({A})=nn 



sinh(Aa - ij) 
\ sinh(Aa - ij + r?) ' 



(3.5) 



For the computation of these expectation values, one needs then to act successively on 
the left with all the elements T^^^^^i {(,m) of the monodromy matrix. One thus has to use the 



expressions of the action of the operators A, B, C, D on an arbitrary state = IljLi C'(Aj). 
In the case of A and D, they are given by |16]: 



N 



N 



N+1 



n sinh(Aj -Xb + v) n+i 
{0\llCiX,)AiXM+i) = aiX^y^^ (0| n C{X 

b=l 



N+1 

n sinh(Aj - Aft) 

3^6 



(3.6) 



N 



N 



(0|nC(A,)^(Aiv+i) 



N+1 n sinh(Aa - Xj + 7?) 



(3.7) 



a=l 



n sinh(Aa - Xj] 



The action of the operator B is more complicated, and it is similar to the successive action of 
A and D: 



N 



N 



N+1 n sinh(Aa - Xk + rf) 



\{C{Xj)B{XN+l) = Y.d{Xa) 



k=l 



a=l 



N+1 

n sinh(Aa - Xk) 



a{Xa' 



N+1 

^ ^h(A7V+l -Xa'+T]) . , - 

n smh(A,-Aa 



N+1 

n sinh(Aj - Xa' + r]) 

(oi n c{x,] 



(3.8) 



Finally, the action of the operator C is free. Recall once more that in the formulas (p.6|)- 
( |3.8| ) the parameters {A} are arbitrary complex numbers (which are not necessarily solutions 
of Bethe equations). In the above sums, the terms containing a(AAr+i) or d{X]y+i) are usually 
called direct terms, while the others are called indirect terms. 
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The action of an arbitrary monomial ^^^^^'^(^i) • • • 7'em,ej„('^m) on the state {ipl can be obtained 
by applying recursively these formulas. This leads to a linear combination of states, 



N 



{0\l[C{Xk)-T,^,4^,)...T,„^,,,jU) = Y.''i{0\ n Cifik), 

k=l i€l k€Ki 



(3.9) 



with some (computable) coefficients a^. Here, sums and products are taken over (multiple) 
sets / and Ki, where the Ki, i G I, are subsets of 1, . . . ,m + N with {fii, . . . , fim+N) = 
(Ai, . . . , Aat, ^1, . . . , S,m)- 



Finally, to evaluate the expectation value ( |3.1[ ), it remains to compute scalar products of 
the type 

TV N 

where Y\f=i B {Xj)\0) is an eigenstate of the transfer matrix, while the parameters {/^j}i<j<Ar 
are arbitrary. The result for (|3.10|) is given by |]l^, (see |]l4| for another proof): 



Proposition 3.1. \1S{ , \2% The scalar product of a Bethe state with an arbitrary state of 
the form l \2.Si ) can be expressed in the following way: 



N 



N 



N 

n sinh(A6 - Ha + r]) 

a,b=l 

Iv 

n sinh(/ia - fib) sinh(Afc - Xa 

a>b 



-det^^''(M[{A}), (3.11) 



/or {Aj}i<j<7v solution of the Bethe equation, and for any set of complex parameters {fij}i<j<N- 
The N X N matrix ^'{{fi}\{X}) is defined by 



^k{{f^}\W)=tiXj,iik)-d{fik)tifik,Xj)l{ 



sinh(/ifc - Xg + r]) 
^jL sinh(^A,. - Xa-T])' 



(3.12) 



Here and further det^y denotes the determinant of an x matrix. Setting {//} = {A} in 
( 3.1C| ), one obtains the square of the norm of the corresponding eigenstate [23|: 



llCiX,)llBiX,)\0)=sinh% n i r\^'^ det^-^^({M)- 



.^ii sinh(Aa - Afe) 



where 



ik 



d{Xj 



N 



+ K{Xj-Xk), 



(3.13) 



(3.14) 



a=l 

with K{X) given in ( 2.15| ). Note that both matrices and can be written in the form of 
Jacobians |14, We also would like to point out that the entries of the matrix are linear 
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combinations of t{\j,^k) and t{l^k^^j)- In the next section we also deal with determinants of 
matrices possessing similar structure. In fact all these determinants are various deformations of 
det t{\j^^k) describing the partition function of the six- vertex model with domain wall boundary 



conditions ||23|. An explanation of this deformation was given in [14| 



Using the formulas (^.2|)"( 3^T^ ), one can compute the normalized average value ( |3.lD on 



the finite lattice. The remaining step is to proceed to the thermodynamic limit. Observe that 
the actions ( |3.6D -(^^) produce sums with respect to parameters {A}. The successive action 

of several E-^' ^ (i.e. successive action of the entries of the monodromy matrix) gives multiple 
sums. Let the parameters \j describe the ground state; then, in the thermodynamic limit, each 
of these sums turns into the integral 

/ fWPtot{\)d\ (3.15) 

{A} ■'^ 



where ptotW is the solution of Lieb equation ( 2.14| ). It was shown also in [15| that in the thermo- 



dynamic limit the ratio of the determinants of ^' and <!>' can be evaluated as a determinant of in- 
homogeneous densities ( 2.16| ). Namely, if in ( ^.11 ) we have {fi} = {^i, . . . , ^n}U{An+i, • • • , Aat}, 



fl{Mptot{Xa)r'det^p{X„^k). (3.16) 



detAr<I>' 

a=l 



In 1 15], it has also been shown that, when acting with the monodromy matrix elements, one 
can get rid of all direct type terms in the thermodynamic limit: the procedure is just to shift 
properly the integration contour of the corresponding A variables. In this way, the expectation 
value of (.^i) . . . T^^^^i^ (^m) in the ground state is represented as a multiple integral in which 
the number of integrals coincides with the number of operators in the product. This quantity 
is called an elementary block. 

This method gives generic answers for the computation of the expectation values of the 
monomials of the type T^jy^iCi) ■ ■ ■ Te^^e'^i^m) corresponding to elementary blocks of the cor- 
relation functions ]l4| , [T5[ | . The problem of the evaluation of the spin-spin correlation functions 
is more involved. Let us consider, for example, the correlation function {erf a^_^_j^) . From the 
solution of the inverse scattering problem (Theorem p.l]) , one obtains the identity 



m m+1 



a; a^^, m ^ (V'l C(ei) • + ^)(^'^) • • n + IV'). (3.17) 



a=2 6=1 

Here {ip) is an eigenstate of the transfer matrix (A + D), and therefore it is straightforward 
to act on the right with nS^^(^ + ^)~^i^b)- However, the determination of the action of 
n^2(^ + is niore involved: clearly after acting with C(^i) on (V'l (or equivalently with 
i?(^m+i) on IV'))) one obtains a sum of states which are no longer Bethe states; therefore the 
multiple action of (A + D) on these states is not simple. In fact, in the framework of the above 
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approach, the product riaL2(^+-^)('^a) would be computed as a sum of 2™~^ monomials, which 
eventually leads us to the sum of 2"*"^ elementary blocks. As already discussed in Introduction, 
this form of the result is not suitable, in particular at large distance m. Therefore, to obtain 
manageable (re-summed) expressions for spin-spin correlation functions, it is essential to obtain 
an alternative and compact evaluation of the multiple action of the transfer matrix on arbitrary 
states. This is the subject of the next section. 



4 The multiple action of the transfer matrix {A + D){x) 

The new results obtained in this section play a central role in the computation of spin-spin 
correlation functions at distance m: we evaluate here the multiple action of the transfer matrix 
^ + D on an arbitrary state, which enables us to solve the problem mentioned at the end of the 
previous section. 

More precisely, let us consider the product 

m 

\{{A + ePD){xa), (4.1) 

a=l 

where xi, . . . ,Xm and [3 are arbitrary complex. Our aim is to find a compact formula for the 
action of this operator on a state (V'l = (0| Ilj^i C'(Aj), where parameters \j are also arbitrary 
numbers. For simplicity, we first consider the case m < N. 

Due to ( |3.6[ ), (^), the action of single operator [A + e^D^ (x) on the state {ip\ can be 
written in the form: 

N N N N 

(0| J] C{X,) (a + e^D) (x) = A(0| J] C(A,) + Ak{0\ J] C{X,) ■ C{x). (4.2) 



j=l j=l k=l J=i 



Here 



and 



N N 

A = a{x) II f{Xj,x) + e^d{x) J] f{x, A,), (4.3) 

3=1 J=l 

smh(x — Afc) smh(Afc — x) 

Recall that a(A) = 1; however, up to the end of this section, we do not use this property, nor 
the explicit form of d{X). All our derivations are based on the use of the equations ( |3.6D , ( |3.7D , 
which in turn are direct corollary of the intertwining relation ( |2.7] ). Thus, the result ( |4.9| ), ( [4.10 ) 
below for the action of the operator (4.1) on an arbitrary state is valid for any model with an 
i?-matrix of the form ( p.ip . 

Let us denote the first and the second type of action in the r.h.s. of ([4.20 by 'direct' and 



'indirect' actions respectively. The corresponding coefficients (O) and (4^) are called direct 
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and indirect terms. The direct action of (^A + e^D^ (x) preserves the state as if this state 
was an eigenstate of the operator. The indirect action leads to the remaining terms, where the 
resulting states differ from the original one: in these states, the argument x of the operator 
(^A + e^Z?) (x) becomes the argument of one of the operators C. 

One can similarly define direct and indirect action for the operator (4J). Namely, the direct 
action of this operator corresponds to the term where the resulting state coincides with the 
original one, which gives 



n + e^D 



a=l 



(dir.) 



rn ^ N AT X 

e''d{Xa)\lf{XaA3)\ 
a=l ^ 7 = 1 7=1 ^ 



(4.5) 



All the remaining terms result from the indirect action. For m < N, a particular kind of 
indirect terms corresponds to the case where the whole set of arguments {x} of the opera- 
tors (j4 + e^-D) (xa) enters the resulting states. The corresponding action is called 'completely 
indirect'. It can be written in the form 



m 



a=l 



(c.—ind.) 



Sm{{x}\{XaJ\{Xa^) {0\l[C{Xa) J] C{Xb) 



{A} = {Aq , }U{Aq_} 



\a\ I — m 



a=l 



(4.6) 



Here the set of the parameters {A} is divided into two subsets: {A} = {Xa+} U {Aq,_}. The 
subset {\a+} is replaced with {x} in the resulting states, and therefore the number of elements 
in this subset is equal m (we have indicated this fact in ( ^ ) by |a+| = m). The subset {Aq_} 
remains in the resulting states. The sum in ( [4. 61 ) is taken with respect to all such partitions of the 
parameters {A}. Following the tradition, we call the corresponding factor S'm({a;}|{Aa4, }|{Aa_ }) 
the highest coefficient. 



Lemma 4.1. The highest coefficient Sm in (4^) can be expressed in the following form: 



5„({x}|{A,J|{A«_}) 



n n sinh(xa - Xb + rj) 

6ea+ a=l 

m 

n sinh(Aa - Afo) H sinh(xa - 

a<b 



Xb) 



det M. 

fc — 1, . . . ,m 



jk 



a>b 



(4.7) 



where the m x m matrix M is given by 



sinh(Aj — Xb + rj) 



Mjk = a{Xj)t{xk,Xj) n f{Xa,Xj)-e^d{Xj)tiXj,Xk) J] /(^j'-^a)]! ■ _ 

sinii(/\'j Xh 71) 

a£a- aGa_ b=l ^ ■' 01/ 
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The proof is given in Appendix 



Thus, we have described two Umiting cases in the action of the operator (4.1): none of the 
parameters {x} enter the final state, or all of them enter the final state. The complete action of 
the operator ( [4.1D contains also all the intermediate cases, when only a subset of the parameters 
{x} enters the final state. 



Proposition 4.1. Let p = m.in{m, N). The action ofYY,^=i + e^D) (xa) on a general state 
written as 



{0\Ylj^iC{Xj), for any sets of complex parameters {xa}i<a<m and {Xj}i<j<N , can be 

m 



a=l 



^ Rni{x,J\{x,_}\{XaM>^ajm J] C(^a) H C{Xb). (4.9) 

n=0 {A}={Aa^}u{Aa_} a67+ bea- 

{x}^{x~f^ }U{x^_ } 
a+l = l7+l=ri 



The coefficient Rn in (4j^) is given by 



Rni{x^+}\{x^.}\{Xa+}\{Xa.}) = Sn{{x^+}\{Xa+}\{Xa_}) 
X n {"(^'^) n f(^b,Xa) n f{h,Xa) + e^d{Xa) f{Xa,Xb) J] /(Xa,A,)|. (4.10) 

Proof. Let us start with the case p = m. Note first that, in this formula, we deal with 
two partitions of the sets {A} and {x}: the subset {A^^} is replaced with the parameters {x^_|_} 
in the final state, and the number of elements in these subsets is equal to n; the remaining 
variables {x^_} do not enter the final state. Once again, the sum is taken with respect to all 
possible partitions of such type. The limiting cases n = and n = m correspond to the direct 
and completely indirect actions respectively. 

One can now use standard considerations of algebraic Bethe ansatz. Namely, let us fix a 
certain partition in ( |4.9[ ) and try to obtain only the state (0| nae7+ ^(^") Hfeea. C{Xb), applying 
the operator (4.1) to the original state (■01- Due to their commutativity, one can re-order the 



operators in the product (4.1) as 



n + e^n) (xa) = n + ^^^) (^-) • n + ^^^) (^«)- (^-i^) 

a=l a£7+ ae7- 

Then it is easy to see that the action of the first group of operators must be completely indirect, 
otherwise one of the parameters rc^, a G 74., is missing in the set of arguments of the final state. 
Moreover, the subset of parameters replaced in the original state {ipl must be exactly equal to 
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{Aq^}: otherwise, at least one of the elements of {Aq,,^} or of {x^_} belongs to the final state. 
Therefore, the action of the first product of operators contributes as 

{^\ n {A + e^D^ (xa) n {A + e^D^ (x„) 

Sn{{x.,+ }\{Xa+}\{>^aJ) (0| J] C{Xa) J] ^(^b) H {A + (x„). (4.12) 

067+ bea- ae7_ 



In its turn, the action of the remaining group of operators (A + e^-D) (xa) in (4.12) must be 



direct, otherwise one of the elements of should appear in the final state. Using (4.5) we 

immediately arrive at ( 4.10| ). 



This proof obviously extends immediately to the case m > N. We just have in this case 
additional direct type actions, and hence no completely indirect term. Thus, the action of the 
operator (|4.1|) on the state is given by (|]|) with the coefficients Rn defined in (p^ . □ 

For our purposes, we will need to specify the arguments {x} of the operator (p|). In the 
inhomogeneous model, we should set them equal to the inhomogeneity parameters: xj = ^j. 
Then, since d{^j) = 0, the part of Rn corresponding to the direct action simplifies (the operator 
D{^j) has no direct action). Thus, the equation ( ^ ) takes the form {p = min(m, A^)): 

m p 

(^1 n + e^D) = J2 5„({e7j|{A.J|{A,_}) 

a=l n=0 {A}={Aa+}u{Aa_} 

{€} = {«7+}U{e7-} 
|cK-|_ I — I'y-)- \ —n 

X n {«(^«) n f^^b,u) n f(^b,^a)\ m n ^(ea) n ^(^b). (4.13) 

Finally, one can also consider the particular case when (-^l is an eigenstate of the transfer 
matrix. Since the parameters {A} satisfy the system of Bethe equations (|2.9| ), one obtains: 



Corollary 4.1. If is an eigenstate of the transfer matrix, that is if the parameters {A} 
satisfy the system of Bethe equations the multiple action of {A + e^D){^j) is given by 



p 



n + e^^) (Ca) = E E ^"({e^JKAaJ) 

a=l n=0 {A}={Aa^}u{Aa_} 

{«} = {«7+}U{e7-} 
|a!^| = |7^|=n 



n n f(^'>,ca) n n /(^^'^«) n n /(^-^^)(oi n ^(^«) n ci>^b).{^.u) 
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Here p = mm{m, N), and 



n n 



n n sinh(ea - A^, + r?) 



Sni^l,...U^l...,Xn)= „ n ^let % , (4.15) 



n 

n sinh(Aa - Xb) U sinh(^a - ^b) 

a>b a<b 



where the n x n matrix Mj^ is 



' sinh(Aj - Aa + 7?) smh(^a - \j + rj) 

Here we have set a(A) = 1. The remarkable property of Sn is that Sn = 6n,o at /? = (see 



Appendix^. Thus, in this case, all the terms in (4.14) with n>l vanish. The subsets {Aq^} 
and {^^_|_ } are empty, and we arrive at 

m m N 

n + ^) = n n f^^^, ^am, (4.1?) 

a=l a=l 6=1 

as it should be for the eigenstates of the transfer matrix. 

5 Generating functional for the correlation function 

In the present and in the next sections, we use the above results to compute the spin-spin 
correlation functions {o'fcr^_^_i). As we are in the framework of algebraic Bethe ansatz, we 
begin with the evaluation of normalized expectation values on the finite lattice with respect to 
an arbitrary eig enstate IV') = nf=i B{Xj)\0). In this case, the relationship between m and N is 
not fixed. However, for the ground state in the thermodynamic limit, tends to infinity while 
m remains finite. Therefore it is clear that we can consider the case m < N only. 

One of the simplest applications of the formulas obtained in the previous section is the 
computation of the generating functional for the correlation function of the third components 
of spins {ala^j^i). Following [|^, IC], we define an operator Qi^m as 

^ m 

Ql,m = ^Y.^l-<Tl). (5.1) 

k=l 

The generating functional is equal to the expectation value (exp(/3(5i,m)), where /3 is some 
complex number. In particular, taking the second derivative of this quantity with respect to (3 
at /? = 0, we obtain (Qf ^). Taking then the second lattice derivative we have: 

^((1 - - af,+i)) = {Ql„,+^) - (Q?,„) - (Qi,„+i) + {QlJ. (5.2) 

Thus, the two-point correlation function of the third components of the local spins can be easily 
extracted from the expectation value of the operator exp(/3Qi^m) in the homogeneous limit. 
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Proposition 5.1. The correlation function (cfo'm+i) in the homogeneous model is given by 



CTi cr, 



l"m+l/ 



2Vi 



(5.3) 



(3=0 



where the symbols and mean the first and the second lattice derivative respectively: 



Vmf{m) = fim + 1) - f{m), Vif{m) = f{m + 1) + /(m - 1) - 2/(m). 



(5.4) 



The expression of the ground state expectation value of the operator exp(/3Qi,m) in the 
inhomogeneous case is 

n=0 ^ ^ p 7=1 b=la=l-'^ 



Here 



and 



xTy„({A},{z}) -detn M,fc({A}|{z}) 



det^ 



/o(Aj,2;fc) 



W^n({A},{z}) = nn 



a=l b=l 



sinh(Aa - Zb + Vj) smh(zfc - Xg + rf) 
sinh(Aa - Xh + v) smh{za - z^ + rf)'' 



M,fc(|A||{.|) - t(z,, A,) + e t(A„ zfc) |i ,i^h(A, - A, + r?) sinh(z, - A, + ^) ' 



(5.5) 
(5.6) 

(5.7) 



T/ie contour C, as in the Lieb equation l \2.14 ), depends on the regime and on the value of 
the magnetic field, and T surrounds only the singularities at the inhomogeneity parameters ^k, 
1 < k < m. 

In the homogeneous limit = r//2, the generating functional of the two-point function can 
he written as 



1 



n=0 



-j-j- C% 

r J=i 



(exp(/3Qi,„)) - J 11 2^i J^- -11 Vsinh(z„ - 2)sinh(A„ + 2 



-Q ( sinh(za + 2) smh(Aa - 2) \ ™ 



xVF„({A}, {z}) • det„ [%({A}|{4) 
where T surrounds the point ri/2. 



(5.8) 



Proof. The operator exp(/3Qi^m,) has a very simple representation in terms of the entries of 
the monodromy matrix. Due to (|3.2| ) we have 



exp(/3Qi,„) = n (a + e^I?) (^a) • H + ^)"' ^^b)- 



a=l 



(5.9) 
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Hence, we can directly apply the equation ( 4.14 ) for the computation of (exp(/3(5i^m)). We begin 
with the derivation of the normalized expectation value of exp(/3Qi^m) on the finite lattice with 
respect to an arbitrary eigenstate of the transfer matrix, = HjLi 

m m 

iM n (A^ pPr)\ (PA n ( a . 

(V'|exp(/3Qi,^)|V^) 



a=l h=l 



The action of the operators {A + Z?)~^(,^b) to the right is trivial: 

m m N 



(5.10) 



+ D)-\i,m = n n r\K,ii>m- 



(5.11) 



b=l 6=1 a=l 

The action of the operators i^A + e^Z?) {^a) to the left is given by (|1|)-(|1|). We obtain 

m N m 



|exp(/3Qi,m)|-^) 



6=1 a=l n=0 {A}={AQ_,_}u{Aa_} 

{C} = {C7+}U{C7_} 
|a+| = |7+l=i 



X n n /(^'"^") n n f^^^^^-) n n /(^-^^) 

aG7_ 6G7+ aG7- feGo_ aGo- 6ga+ 

AT 

(oi n c{ia) n c^(A6) n ^(a,)io) 
(oi n c{\,) n Bixm 

i=i i=i 

In fact, for the remaining part of the calculations, we can use the formulas of Section ^. Using 
( 3.11| )-( pj.l4| ) we find the ratio of the scalar products in the r.h.s. of ( ^.121 ) and finally arrive at 

'""t^r""' ^ s E n n n n 

n=0 {A}={Aa+}u{AQ_} b=laGa+ ae7_ feg^,,. 

{«} = {«7+}U{«7_} 
|a+| = |7+|=n 

n n sinh(ea - Afe + 7?) sinh(Ab - + f?) ^ ^ 

"^7+tea+ detjv^^({g^^}U{A»_}|{A}) 

n sinh(Ca-6) n sinh(A,-A, + r?) ■ A ' det,v^'({A}) ' ^ 

<>>i'e7+ a,bGa+ *^^+ 

The representation ( 5.13 ) gives us the expectation value of the operator exp(/3Qi^m) on the 
finite lattice. 

Up to this stage, our derivation was purely algebraic. Now we should proceed to the ther- 
modynamic limit. In spite of the fact that this limit strongly depends on the phase of the 
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model, one can present the final result in a quite general form. First, using ( |3.16| ), we have in 
the thermodynamic limit 

det^^'({e^Ju{A„_}|{A}) -p. n ^\-iA , a ^ (r,,A^ 

= _[_[ (Mptot(Aa)) det„p(Aj,^fc), (5.14) 



det7v$'({A}) 



where \j G {Xa+} and £ {^7+}- Second, we should get rid of the sum with respect to 
partitions in ( ^.131 ), replacing them with integrals. 

Consider first the partitions of the set {^}. Let J-{zi, . . . , Zn) be a symmetric function of n 
variables Zfc, analytical with respect to each argument in the vicinities of {^}. Then 

{«}={«7_}u{e7+} ae7- b67+ 

7_|_ I =n 

n n 

n n sinh(za - Zb) 

=;^/n|^nn/(--^^)W^^^ ^(w). (5.15) 

r i=i »=i ''=1 n n smh(za - Zfc + r/) 

a=l b=l 

Here the contour T surrounds the points ^^i , • • • , and does not contain any other singularities 
of the integrand. In particular, when all — > r//2 (homogeneous limit), one can chose V as an 
enough small circle around ri/2. Observe also that one can easily take the homogeneous limit 
in the r.h.s. of ( 5.15| ), while the existence of such limit in the l.h.s. is not so obvious. Thus, 



we can replace the sum with respect to the partitions of the set {^} with the set of contour 
integrals. In fact, we could do this already on the finite lattice. 

As for the sum of the partitions of the set {A}, we essentially use the properties of the 
thermodynamic limit. Let now J^(Ai, . . . , A^) be a symmetric function of n variables, vanishing 
at Xj = Xk, j, k = 1, . . . ,n. Then, 

1 ^ N N 

Yl •^({Aa+ }) = ^^J^ X] ■ ■ ■ X] -^(^ii ' • • • ' ^in) 

{A} = {A„^}U{Aa_} ■ ji = l j„ = l 



\a_i_ \ —n 



nl Jc Jc 



dXiptot{Xi) ■■■ I dXnPtot{Xn)T{Xi, Xn), as M ^ oo, (5.16) 



when C is the contour involved in the Lieb equation ( p. 14 ) 



Thus, the sum with respect to partitions in the thermodynamic limit transforms into 2n 
integrals, and we obtain the expectation value ( |5.5| ) of the operator exp(/3(5i,m) in the ground 
state. □ 

We would like to mention that, as for the elementary blocks of correlation functions ob- 
tained in ||l5||, the integrand in (|5.^) consists of two parts: thermodynamic and algebraic. The 
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thermodynamic part depends on the phase of the model (A and h), but not on the considered 
set of local operators; it includes the determinant of densities and the integration contour C 
for the variables {A}. On the contrary, the algebraic part, which includes the remaining factors 
of the integrand, does not depend on the phase of the model, but only on the particular set of 
operators the correlation functions of which one wants to compute. 



Nevertheless, there exists a principal difference between the representation (5^) and the 
multiple integrals for the elementary blocks: 

Each of the elementary blocks involved in the construction of JlaLi + e^-D) (^a) contains 
exactly m integrals. However, some of these integrals, which correspond to the action of -D(^a), 
are taken over the contours C (see ( p. 14]) ), while the other integration contours (corresponding to 
A type actions) are shifted. The reason of this difference is that, due to the fact that d{S,a) = 0, 
the operator -D(^a) has no direct action, while A{^a) does have. The shift of the contours allows 



us to get rid of the terms produced by the direct action of A{S^a) (see [15| for more details). 

On the contrary, in the equation (^.^), the integrals with respect to all {A} are taken over 
the same contour C. Moreover, the number of these integrals is not fixed, but varies from to 
m: the term n = m in (|5.8| ) corresponds to completely indirect actions of A and D; the terms 
with n < m contain direct actions of A, and in particular the term n = describes the direct 
action of the whole product A{^i) . . . A{^m)- 

From this observation, one can easily understand how to obtain the representation ( |5.8D 
from the formulas for the elementary blocks. First of all, one should move back all the shifted 
contours in each of the 2™ constituent elementary blocks of the product HoLi + e^-D) {^a)- 
Hereby one crosses the poles of the densities p{\, ^), which produces the terms where the number 
of integrals is less than m. Then, one needs to gather the terms with the same number of 
integrals and symmetrize all the obtained integrands with respect to {A}. The symmetrization 
produces the sum over the partitions of inhomogeneities {^}, which can be effectively taken 
into account by the set of auxiliary z-integrals. Generically, this way meets extremely serious 
practical difficulties. However, for some simple particular cases, it can be successfully applied. 
In Appendix we illustrate this method by considering the limit /3 — > cx3 in the equation (|5.8D . 
Then we have 

m m 

hm e-^™ exp(/3Qi,^) = TT D{ia) \\{A + D)-\i,). (5.17) 

a=l 6=1 

Hence, we obtain the correlation function corresponding to the emptiness formation probability 



which was considered in Section 4 of |15|. On the other hand, in the equation (|5.8| ), only the 



term with n = m survives, and the determinant det simplifies: 

hm e-^"'detmM,k= f\ ^^^M^'' " + ^) det^[t(A„ z^)]. (5.18) 
is^co \^ smh(2;a - h + V) 

In Appendix |^, we obtain this result by symmetrization of the corresponding multiple integral 
obtained in |jl5|. Finally we observe that, due to the fact that det„Mjfc = 6no at /3 = 0, the 
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expectation value of the identity operator is equal to 1. This obvious fact would be highly 
nontrivial to prove from the multiple integral representations for elementary blocks. 



6 Spin-spin correlation functions 

In the previous section, we have obtained a multiple integral representation of the correlation 
function {(jfa^^-^) from the generating functional (exp(/3Qi^m))- It has been computed from 
the average value in the ground state of the product of m commuting operators {A + e^D){^a) 
for a = l,...,m. However, for general n-spin correlation functions, we have to deal with 
ground state average values of products of non-commuting operators. Indeed, a generic A;-point 
correlation function can be written as 

(V'ln^-ri^)' (6-1) 

where mi < m2 < . . . < rrik is an ordered set of k sites on the lattice. Using the solution of 
the quantum inverse scattering problem, it can be reduced to the evaluation of the following 
average value: 

m2 — 1 7713 — 1 

ii^K^A^^m,)- n + n {A + D){Ca)...T,„,'^{Um, (6.2) 

a=mi+l a=m2+l 

To compute such a correlation function, the strategy is the following. As usual, we act with 
the operators to the left. The action of the shift operators (products of transfer matrices) is 
given as in (|4.9| ). The action of the isolated elements of the monodromy matrix T^-^e'XCmj) for 
j = I, . . . ,k is computed as in the elementary blocks: in particular, in the thermodynamic limit, 
we get rid of all direct type terms by shifting the corresponding integration contours (see |jl5|). 
This helps the final expression to be as compact as possible. We apply below this procedure to 
the case of the spin-spin correlation functions. 

Let us start with (o'lC^+i). 

It is possible to evaluate this quantity by computing directly the expectation value 
(VkiV^+il^) = mA-D)i^i) ■ l[(A + D){Ca) ■ {A-D){U+i) ■ n {A + D)-\i,m. (6.3) 

a=2 b=l 



First, we has to to act with {A — D){^i) on the eigenstate (^/^l using ( p.6D , (3.7); then, we can 



apply ( 4.131 ) for the action of the product 0)^12 (^ + ^)(^a)5 and finally again use ( |3.6| ), ( |3.7D for 



the action of (^4 — D){^mj^i). This method does not meet some new principal obstacles, but the 
final answer has a bit more complicated form. In fact, in this case, one obtains a sum similar 
to (|5.8| ), but in which the n^^ term consists of four summands corresponding to the action of 
the operators A and D in the points and ^m+i- Therefore the evaluation of (ciC^+i) via 
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the generating functional seems to be the most preferable, at least from the stand point of its 
compactness. 

However, in particular cases, other representations for the correlation function (crfo"^^^) 
may present certain advantages, even if these representations have slightly more complicated 
forms. Note that generically the correlation function (cfa^^^) contains the term (u^)^ which, 
due to the translation invariance, does not depend on the distance m. However, at zero magnetic 
field in the massless regime, the magnetization (cj^) is zero and the correlations of the third 
components of spin should be a decreasing function of m. On the other hand, each term of 



the sum (5.8), even after taking the lattice and /3-derivatives, still contains a part which does 
not depend on the distance m. It follows from ( |5.2| ) that eventually these constant terms give 
1/2. For the purposes of the forthcoming asymptotic analysis, it would be desirable to obtain 
a representation of (cfCm+i) without any constant contribution at zero magnetic field. Such a 
representation, of course, is provided by (|6.3D, but we can also consider a certain modification 
of the generating functional (exp(/3Qi^m))i namely {exp{(3Qi^m)crm+i) ■ Indeed, it is easy to see 
that 

1 d r 

' exp(/3Qi,m) - exp(/?(52,r 



0-1 (T^ 



and, since the magnetization is zero at h 



d 



f3=0 

in the massless regime, we obtain 

i) , h = 0, |A| < 1. 



(6.4) 



(6.5) 



/3=0 



Proposition 6.1. The ground state expectation value (exp(/?(5i,m)f'"m+i/ 
case is given as 



for the homogeneous 



(exp(/?Qi,m)o-^+i) = 



-1 



(n!)2 

n=0 ^ ^ P J 



ndzj_ 
97r7 



c 



n 

a=l 



sinh(za + 2) sinh(Aa - 2) 
smh{za - ^) sinh(Aa + f) 



;iy„({A}iM)detJM,,({A}I{z})j J\ ''''^^^'^ ^ 



dX 



n.+l 



n 



sinh(A 



n+l 



C 



L sinh(A„+i - \a 



^}-^ sinh(za - 2) 

, TT sinh(An+i -Zg + rj) 

V) J "'''^ _i sinh(An+i - Xa + v) 



T]) 



det 



n+l 



(6.6) 



p(Aj,Zfc) 



Here one should set Zn+i = ??/2 in the last column of the determinant of the densities. C and T 



are as in Proposition -SA, while the shifted contour C for Xn+i is such that CU (— C) surrounds 
the points {z}, in which the functions p{Xn+i,Zk) have simple poles, but it does not contain 
other singularities. For zero magnetic field one has 



(exp(/3Qi,m)a^+i) = ^ 



1 



n=0 



(n!)2 



ndzj^ 
9.Tri 



a=l 



sinh(2;a + 2) sinh(Aa - ^) 
sinh(za - 2) sinh(Aa + ^) 
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n( sinh(Aa - V. 



\smli{za 

a=l ^ 



n 



smh(A„+i - Za) 



^J-^ smh(A„+i - Xa) ^J-^ sinh(An+i - Xa + r]) 



xdet„ 



detn 



+1 



p{Xj,Zi), . . . ,p{Xj,Zn),p{Xj, 2) 



(6.7) 



Proof. We ommit parts of the proof which coincide with parts of the proof of Proposition 



5.1| . Instead, we focus our attention on certain pecuUarities. 

The normaUzed expectation value of exp(/3Qi^m)o'm+i on the finite lattice is given by 



exp{(3Qi^rn)cr: 



m+ll 



m m+1 

\W{A + ePD){ia)-{A-D){U+i)- n {A + D)-\i,m 

a=l 6=1 



(6.8) 



m+l 



After the calculation of the actions of the product W + D) ""^(^b) on the state and of 

b=i 

\Y^=i{A + e^D){^a) on the state {ip\, one has to act with {A — D){S^rn+i) on the resulting states 

{ma^,^C{^.)Utea^C{X,). 

The action of D{S,m+i) gives 

(oi n ^(^'^) n ^(^^) • ^(^-+1) - ^ 



E 



sinhfA/ — fm+i. 



n /(^^'^-) n /(^-^^) 



(6.9) 



Here we have used that {A} satisfy the system of Bethe equations. Then, one has to compute the 
normalized scalar products of the obtained states with \^|J) and proceed to the thermodynamic 
limit. We obtain: 

77T, 71, Th TTX 

ie..m..^mu..)} = E / n Si X ■ n n 

n=0 ^ ^ p j=l -"^ b=la=l''^ "'^"'^ 



:14^„({A}|{.})det jM,,({A}|{z})j J] 

L J ^_-|^ smni^Zfl sm+i / 



X / dXn+l ■ Jl 



sinh(A„+i - Za + ri) 



■ detr 



p{Xj,Zi),.. . ,p{Xj,Zn),p{Xj,Cm+l) .(6.10) 



ic sinh(A„+i - Xa + rj) ""^^ 

Observe that, in comparison with ( |5.8| ), the integrand in ( |6.10 ) contains additional factors. 



Moreover, we see that the sum over £ in (|6.9| ) produces one more integral with respect to A^+i 
in the thermodynamic limit. 
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The action of the operator A{(^rn+i) is more compHcated. First, it contains terms similar to 

^ smh(r!i'- A.) n /(^-^^) n /(Aa,A,).(0|C(W) U C{i.) H C{X,). 

In the thermodynamic hmit, the contribution of these terms turns into integrals of the type 
( 6.1C| ), where in the last line one should make the replacement 

-pr sinh(An+i - Zg + r]) ^ _ fV sinh(An+i - Zg - rj) 

ijL sinh(An+i - Xg + v) sinh(A„+i - Xg-rj)' 

However, the action of the operator A((^m+i) contains also the direct term 

n n /(Aa,w)-(oi n ^(^-^ n ^(a^-)' 

ag7+ aeo- ae7+ bea- 

and terms where the argument £,m+i exchanges with one of {^-y^}: 

^ sinhr!.'- 6) n /(?-^^) n /(Aa,e.)-(oic(e„+i) n c{^a) n c{x,). 

These terms give also non-vanishing contributions to the expectation value (exp(/?(5i,m)o'm+i)- 
However, in the thermodynamic limit, one can withdraw them by shifting the integration con- 
tour for A„+i in complete analogy with the method used in for the elementary blocks. This 
gives us ( |6.6D . This representation is valid in arbitrary regime of the model. For zero magnetic 
field, the original integration contour C is the real axis for | A| < 1 and the interval [^7r/2, —111/2] 
for A > 1. In both cases we can choose C = C + r] (recall that in the massless regime r] = —iC, 
C > 0). Changing then A„+i with A„+i + i], we arrive at (|6.7|). □ 



The representation (|6.7| ), as it was expected, has a slightly more complicated form than 
(|5.8| ). However, after the lattice and /J-derivations, the sum (16.71) does not contain any constant 
contribution for |A| < 1. This fact may play an important role for the asymptotic analysis of 
the correlation function (o'fo'm+i)- 

The remaining two-point functions are ^^'^ i^t^m+i)- ^n the finite lattice these 

two quantities are given by 



m+l 



|CT1 0-^+1 1^) _ a=2 6=1 



m m+l 
\(T^ (T^+l\lp) _ g=2 b=l 
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(6.11) 



(6.12) 



It is clear that our method can be apphed without significant changes for the calculation of 
( 6. lip , ( |6.12| ) as well. Therefore we present here only the final result. Moreover, for simplicity, 
we first consider the case of zero magnetic field (for both massive and massless phases), when 
the correlation functions (ffo-m+i) ('^I'^m+i) coincide. 

Proposition 6.2. The ground state expectation value for the homogeneous case at 

zero magnetic field can be expressed as 



m— 1 



1 



n=0 



nUn + 1)! 



r J=i 



d 



n+2 



C 



n+1 
a=l 



sinh(z„ + 2] 
sinh(2;a - 2] 



sinh(An 



+1 



An 



+2J 



n 

a=l 



/n+1 

n sinh(A„+i - Za + ri) sinh(A„+2 - ^aj 

a=l 

n 

n sinh(A„+i - Aa + r?) sinh(A„+2 - K) 

\a=\ 



sinh(Aa - W 
sinh(Aa + \) 



M^n({A},{4) 



xdet„+iMjfc ■ det„+2 [p(Aj, zi), . . . , p(Aj, z^+i), p(Aj , 2)] , (6.13) 



where the contours C and T are defined as in Proposition \5. 1\ . Here the analog of the function 

Wn{{X},M) 

n n+1 

n n sinh(Aa - Zb + r]) sinh(z6 - Xa + r]) 

^"(W' = ' (6-14) 

n n smh(Aa - Afo + r?) n n smh{za - Zb + ij) 

a=l 6=1 a=l b=l 

and the (n+1) x (n+1) matrix M has the entries 

tCj .( \ \ .(\ A TT sinh(Aa - \j + r?) sinh(Aj -Zb + ri) 

M,, = t(z„ A,) - t(A„ ..) n n ,i,h(z.-A,+,) ' ' - 

and Mn+i,k = t{zk, ^) for j = n + l. 

Proof. The only difference between ( |6.12D and the expectation values considered above is 
that now we deal with the operators C and B, and thus use ( |3.8D for the action of B (recall 
that the action of C is free). The rest of the computations is mostly the same. In analogy with 
(|6.7D , the integrals with respect to An+i and A„_|_2 describe the action of the operator i?(^m+i) 
in ( [6.121) . The direct action of B{^m+i) is taken into account by the shift of the integration 
contour for A„+2- D 

It is clear that this result can be easily generalized for the case of non-zero magnetic field in 
complete analogy with (^). In particular, for the correlation function {afa^_^_^), one should 
replace the variable Xn+2 in ( 6.13| ) with A„+2 + ri, and choose for this variable the shifted 
integration contour C. 
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Conclusion 



The main result of this paper is a new multiple integral representation for the spin-spin cor- 
relation functions at lattice distance m of the XXZ Heisenberg chain in a magnetic field. In 
particular, it gives generically an effective re-summation of the corresponding 2*" elementary 
blocks as the sum of only m terms, each containing the distance as the power m of some simple 
function. Hence, our method opens the possibility of the asymptotic analysis of the spin-spin 
correlation functions at large distance. It will be shown in a separate publication that it also 
leads in a direct way to the known answers at the free fermion point A = 0. 

It should also be noted that the compact formula for the multiple action of the transfer 
matrix operator, for any values of the spectral parameter and on arbitrary quantum state, is 
central in our result. It can be used to compute multi-spins correlation functions. It contains 
in particular the possibility to act on any quantum state with generic conserved quantities 
responsible for the quantum integrability of the XXZ Heisenberg chain, for example with the 
Hamiltonian itself. Therefore, this formula is also the key to the dynamical correlation functions. 
Note finally that this result depends only on the general structure of the i?-matrix, and thus 
can be generalized to other models admitting quantum inverse scattering problem solution , 
like the integrable Heisenberg higher spin chains [pi] ]. 
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A The highest coefficient 



Let 



n n 



n n sinh(xa - fib + i]) 




6=1 a=l 



• detnMjk, (A.l) 



n 



n sinh(;Ua - /Ub) sinh(xf, - Xa) 

a>b 



where the n x n matrix Mj^ is 

N 



N 



sinh(//j — Xb + r]) 



Mjk = a{fjLj)t{xk,fij) f{fJ.a, fJ-j) - e'^d{fjLj)t{iJ.j,Xk) 




a=n+l 



a- 



,=n+l 6=1 
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First of all we prove an auxiliary lemma, establishing the recursion property of the function 5^- 
Lemma A.l. 

, Xn\^l ) • • • ) /^n I A'n+l i ■ ■ ■ i ^J'N) 

n 

= ^5'„_i(xi, . . . ... ,l2n\tJ'n+l, ■ ■ ■ ,fJ-N,fJ'l) (A.3) 

(=1 

(n-1 N n-1 N \ 

aifJ-l)g{Xn,Hl)Y[fiXa,fJ-l) Y[ fi^J'a,^J'l) + e^d{^J,l)g{nl,Xn)Y[fi^J'hXa) /(W,A^a)|, 

a=l a=n+l a=l a=n+l / 

where the symbol fii means that the corresponding parameter is ommited in the set ^i, . . . 
Proof. Consider an auxiliary contour integral: 
If duj 



-Sn{xi, . . . ,a;„_i,u;|//i, . . . ,;Un|//n+l, • • • (A.4) 



27ri J sinh(3;„ — uj 

The integral is taken with respect to the boundaries of a horizontal strip of the width in. 
For instance, one can take for the lower boundary = wq, and for the upper boundary 

9(0;) = ujQ+iiT. Hereby loq is an arbitrary real number satisfying the conditions ojq 7^ '^{xn+iT^k) 
and cjQ 7^ Q^l^i + «vrA;), where j = l,...,n, k € Z. Obviously, the integrand decreases as 
exp(— 2|a;|) at w ^ ±00. Moreover, the integrand is a periodic function of uj with the period 
iiT. Thus, 1 = and, hence, the sum of the residues inside the contour vanishes. The pole 
at sinh(ti; — Xn) = gives us the term in the l.h.s. of ( [A.3| ). On the other hand, the only 
singularities of the function S„ are simple poles at sinh(a; — Uj) = 0. The residues in these poles 
give us the r.h.s. of ( A.3| ). Thus, the lemma is proved. □ 



Proposition A.l. The highest coefficient of the completely indirect action is equal to the func- 
tion Sn at /ii, . . . ,/i„ = {\a+} and . . . = {Aa_}. 

Proof. One can use the induction with respect to n. For n = 1, the equations ( ^.7] ), ( [4. 81 ) 



give us exactly (4^). Let the highest coefficient have the form (L7), (|4.8| ) for n — 1. Then we 
have 



a=l 



s'„_i(xi,... |{A„J1{A„_}) 



(c.-ind.) {A}={Ac<^}u{Ac<_} 

I Q!_|_ I —n — l 



n-1 



■.{0\l[C{xa) n C{Xb)-[A + efD 

a=l bGa-. 



[Xr, 



(A.5) 



(c.—ind.) 
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Now, acting with the last operator, we need to exchange Xn with one of the A G {Aq }. This 
gives us 



a=l 



5„,_i(xi, . . . ,a:ri-i|{Aa+}|{AQ,_}) 



(c.-ind.) {^}={^a_,_}u{AQ_} iGa_ 



\cti l—n — l 



n-1 



n-1 



a(AO<?(x„,AOn/(^«'^0 n f{>^a,Xi) + e^d{Xi)g{\i,Xn)llf{Xi,x,) J] /(Ai,AJ 



a=l 



a=l 



a=l 66q!_ 

b^l 



For each fixed partition, one can define the new sets 

{a„y} = {a«+}uA/, 

{A„'_} = {A„_}\Az. 



Then, ( |A.6| ) takes the form 

n 



a=l 



(c— ind.) 



(A.6) 



(A.7) 



J2 J2 ^n~l{xi,...,Xn-l\{Xa'^}\M{K'_}^>'l) 



I 



n-1 



n-1 



a(A05(xn,Az)n/(^a'Az) n /(Aa,AO + eMA/MAz,x„) J]/(Az,2;,) J] /(Ai,A„ 

a=l aea'_ a=l aGa'_ 



hea' 



(A.8) 



a=l 



Identifying in ( |A.8|) and (|AJ) Xi = m, {Xa'^] = fii,...,fin and {Aa'_} = . . . , //at, we 

come to the conclusion that the coefficient at the state (0| na=i C'(xa) HftGa' ^(Afe) is exactly 
Sn{xi, . . . ,x„|{Aq' }[{A„' }), what ends the proof. □ 
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B The properties of the function Sj^ 



Proposition B.l. The matrix Mjk{6.i, . . . , ^n|Ai, . . . , A^) ) at (3 = possesses the eigen- 
vector 

n ( n 



(B.l) 



a=l 



\ 



a = l 



with zero eigenvalue. 



Proof. The action of Mjk{ii, . . . , . . . , A„) at /3 = on the vector ( [B.l[) can be written 

in the form 



■,{-) -TT sinh(Aa - \j + rj) sinh(Aj -jg + ri) 



k=l 



L sinh(Aj - Xa + ij) sinh(^a - Xj + ij) ' 



(B.2) 



where 



G 



(±) 



E 



sinh 1] 



n sinh(^fc - Aa) 

a=l 



^ sinh(efc-A,)sinh(e,-A,±r/) ^-^^^^^ _ 



(B.3) 



a = l 



To find gJ^^ we consider a contour integral similar to the integral in Lemma |A.1| 



(±) 



sinhry 



(B.4) 



27ri 7 sinh(a; — Aj) sinh(ti; — Xj ± ??) sinh(a; — ^a) 

Just like in the Lemma A.l , the integral is taken with respect to the boundaries of a horizontal 
strip of the width in. Due to the periodicity of the integrand and its vanishing at w — > ztco, 
we conclude that 1^^^ = 0, and thus that the sum of the residues within the contour vanishes. 

The sum of the residues at sinh(ci; — ^fc) = gives G^f^^ ■ In addition, we have one more pole at 
sinh(t(; — Xj ± r/) = 0. Combining all together we find 



G 



(±) _ ^ sinh(Aa - Aj ± rj) 



(B.5) 



" ^J-^ sinh(^a - Xj ± rj) 

Substituting this into (|B.2| ) we obtain Ylk=i ^jk^^k = 0, and the Proposition is proved. □ 
Thus, at /3 = 0, the determinant of the matrix Mj/. vanishes for n > 1 and, hence, 5„ = 5no- 
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C Symmetrization of an elementary block 

The multiple integral obtained in [|l5| ] for the emptiness formation probability r(m) on the 
inhomogeneous lattice has the form 

n sinh(ea - ib) 

a<b 



C 



where 

m I j—1 m I 

n < n sinh(Aj -Ck + ri) U sinh(Aj - ^k) \ 

/({A}, U}) = '^^^ ^ '-^ L (C.2) 

n sinh(Aa - Afe + r?) 

a>b 

Clearly, due to the factor detm[p(Aj, ^a;)]) the symmetrization of the integrand with respect to 
all {A} is equivalent to the alternating sum of -f({A},{^}) with respect to the permutations 

(T : Ai , . . . , Am ^ Ao-(l) , . . . , Xa{m) ■ 

Proposition C.l. 

W/({A.}, {^D = Zm({A}, {e}), (C.3) 

where 

7 aw sew TT TT s™h(Aa - gb) sinh(Aa - + 7?) detm[^(Aj, g^)] ,s 
^™({A}, m) = ii ii ,inh(A, - A, + r?) . (^-^^ 

Proof. It is convenient to introduce new variables Xj = e^^^ , = e^^-', q = e^. Then ( p.3| ) 
takes the form 

n S n - IVk) n (a^^O) - Vk) \ 

m>a>fe>l 

and 

-1 _ _ 

/ (a^a - yb){q'^xa - qyb)\ 

\a=l / \m>a>6>l / a=l b=l 



m m 



zu{x},{y}) = [llxA\ n (ya-yt)] UU 



{q - qxb) 



xdetr 



yk){q~ 



yk) 



(C.6) 



Observe that ( |C.5|) holds for m = 1. Suppose it is valid for m — 1. Let us consider the properties 
of the both sides of (p.5| ) as functions of i/m- Obviously, these functions are polynomials 
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of Urn of m — 1 degree (the poles in i/m = Va in the r.h.s. disappear due to the zeros of 
the determinant in the same points). The coefficients of these polynomials are antisymmetric 
functions of parameters {x}. Thus, in order to prove ( |C.5D , it is enough to compare the values of 
both sides of this equality in ym = Xa, a = 1, . . . , m. Moreover, due to the antisymmetry of the 
coefficients with respect to {x}, it is sufficient to consider the case ym = Xm- For ym = Xm the 
determinant in ( p.6| ) reduces to the product of the last diagonal element by the corresponding 
minor. Extracting all the dependency on Xm and ym, we obtain for Zm'- 



Zm{{x}, {y}) 



— T ^ Zm-l{{x^Xm},{y^ym})- (C.7) 

[q '-Xm-qxa^ 



Vm — a=l 

m—1 1 



Consider now the l.h.s. of ( C.5 ). Each term of this sum contains the product YYjl=L ixa{j) — Um)- 
For ym = Xm, this product does not vanish if and only if ^^-(m) = Xm- Hence, in this case, we 
need to sum up only with respect to the permutations of the m — 1 variables xi , . . . , Xm-i , while 
Xm remains fixed. Denoting these permutation as a' , and extracting again the dependency on 
Xm and ym, we obtain for the l.h.s. (|C.5|): 



m—1 f i^l m—1 

^ n { UiQ~^x^>(^j) - qyk) n ixa'ij)-yk) 

{Xa-ym){q Xm - qya) i\pf<7')-?=^ l'==l k=j+l 



M {q^^^Xm - qXa) ^ n ^q'^X„,(a) - qXa'ib)) 

" m-l>a>b>l 

(C.8) 

Due to the assumption of the induction, the sum with respect to the permutations of the 
variables xi, . . . ,Xm-i gives Zm-i{{x ^ Xm}, {y 7^ ym})- Then, comparison of ( |C.8| ) and ( |C.7| ) 
completes the proof. □ 



Thus, after symmetrization of the integrand, the multiple integral ( C.l ) for the emptiness 
formation probability takes the form 

rim) = —[ ^^^(Waa) det^[p(A„efc)]rf'"A. (C.9) 
^ ^ nsinh(ea-e6) 

a<b 



ml Jc 



On the other hand, taking the limit /? — > 00 in (^.5|), we obtain 



xWm{{X}, {z}) ■ detm[t{Xj, Zk)] ■ detm p{>^j,Zk) 



(C.IO) 



Taking the contour integrals with respect to all Zj, we immediately arrive at (C.9). 
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